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Abstract 

The ability to pull is one of the primary distinctions between an automobile and a tractor. The tractor slippage 
has a significant impact on the pull. It is critical to keep slippage in the 8 to 15% range to achieve improved 
tractive efficiency, i.e., using more pull with the available power. As a result, slippage has a higher impact on 
the performance of a tractor-implement combination. Currently, slip is measured in the field using manual 
intervention. To solve this issue, a MATLAB Simulink model was developed to forecast the slip for a two-wheel 
drive tractor, and a micro-controller-based wheel slip measurement device was utilized to measure and 
validate the slip. In addition to saving time and costs, predicting slippage would lessen the laborious work 
that humans do during field testing. The developed model could also predict the actual speed of operation and 
the field capacity. The developed model predicted the slip within +4.03 to -4.55%, actual speed of operation 
within +3.83 to -2.70 % and field capacity within +3.44 to -7.09% relative deviation. The coefficient of 
determination was found close to 0.98 for all the three performance factors. The developed model is user 
friendly and efficient. This could help researcher to forecast the performance parameters without going to 
field. 

Keywords: MATLAB Simulink model, Tractor, Implement, Slip, Wheel slip measurement. 


1. Introduction 


Amidst the increase in oil prices and labor expenses, 
considerable attention has been devoted to enhancing 
the productivity of tractor-implement combinations 
during field operations. Slip, as highlighted by 
Wismer and Luth (1972) and Brixius (1987), plays a 
significant role in optimizing tractive performance. 
[l, 2]. As tractive efficiency is optimal within this 
range, maintaining wheel slip within a specified 
optimal range ensures the tractor operates at 
maximum efficiency. [3]. There is no instrumentation 
in the current tractor system to display or measure 
this slip. In field testing, it is done manually. Despite 
the fact that, technology has been developed to 
measure the same. Wherein, two speed sensors are 
fitted onto front and rear wheel of the two-wheel 
drive tractor [4, 5, and 6]. In this study similar 
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methodology was used to measure slip. Since slip is 
a major design factor for tractor-implement 
performance and current tractors doesn’t have the set 
up to measure slip on the go hence field-testing 
methodology is the only way to identify slip. This 
methodology is labour intensive wherein, the 
marking of the field is required and manually 
counting the number of revolutions takes time. It is 
iterative process and different implements are tested 
to arrive at the result. Hence to overcome this 
problem, a MATLAB Simulink model-based design 
was developed that could predict the slip of the 
tractor while in operation with good accuracy. Since 
actual speed of operation and actual field capacity are 
also important parameters and it can be predicted if 
slip is predicted accurately. Therefore, these two 
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parameters were also predicted and included in this 
study. In order to save time and cost and also lessen 
laborious work, this study intent to predict slip, actual 
speed of operation and field capacity which are an 
important factor for performance of the tractor while 
on operation and it was validated by measuring slip 
via instrumentation. 
2. Theoretical Considerations 

To create a MATLAB Simulink Model-based design 
for forecasting slip, actual speed of operation, and 
field capacity, and the following theoretical and 
empirical equations were used. 


2.1 Traction Prediction Models 
Tractive performance of agricultural tractors was 


Table 1 Recommended Traction Equation Coefficients 
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developed using the traction prediction equations of 
Brixius (1987). The following traction equations are 
taken into consideration in this study. 


Mg = AA - ae: x (1-e7435) +A, (1) 
= os Q) 
MRR, = 3+ Ay + sexs (3) 


Where, gis gross thrust force coefficient, B,is 
Brixius number, S is slip (%), b is wheel width (in.), 
dis wheel diameter (m), & is tyre deflection, h is tyre 
section height (m), R,is rear wheel dynamic weight, 
MRR, is Motion resistance ratio and Aj, Ap, A3.....A7 
are coefficient’s of traction. The coefficients are 
shown in Table 1 [1, 2]. 


Traction Coefficient Al A AS At As As AD 
Type of tire 

Bias ply 0.88 0.1 75 0.04 5 3 1 
Radial ply 0.88 0.1 8.50-10.50 0.03-0.035 3 3 0.9 


2.2 Reaction Forces in Tractor 


Force 


—— Distance 


Figure 1 Forces Acting On Tractor-Implement Combination 


Considering forces and moments in Figure 1, the 
dynamic wheel reactions were derived as in equation 
4and 5. 
pw Weep Keg) (Wim epi tHatLtep)-D¥a 4 
ad L-e,+er @ 
Rp = (W, + Wn + Py) — Rr (5) 
Where, CG; is centre of gravity of implement, CG; is 
centre of gravity of tractor. D is draft (kgf), d is depth 
of operation (m), e, is eccentricity for front tire (m), 
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er is eccentricity for rear tire (m), F is gross thrust 
force (kgf), Hg is distance between rear axle center 
and link point of the tractor(m), L is wheel base of the 
tractor (m), Pyis vertical force (kgf), Reis front wheel 
reaction (kgf), Rr is rear wheel reaction (kgf), TF¢ is 
towed force front (kgf), TF; is towed force rear (kgf), 
Wn is weight of implement (kg), W; is weight of the 
tractor (kg), Xcgi is distance between link point and 
center of gravity of implement (m), Xcgt is distance 
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between rear axle center and center of gravity of 
tractor (m) and Yq is center of resistance distance to 
ground surface (m) [4,5,6]. 

2.3 ASABE Draft Prediction Model 
Tn the software calculation, the existing ASABE draft 
model was employed to forecast the draft. Equation 6 
expresses it [7]. 
D= F(A+ BY +CVz)xWxd (6) 
In the given equation, D represents the draft in 
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kilograms, while F_j denotes the dimensionless soil 
texture adjustment, with j taking values of 1 for fine 
soil, 2 for medium soil, and 3 for coarse soil. 
Parameters A, B. and C are specific to the machine 
[8] . W represents the width in meters or the number 
of tines, d stands for the depth of operation in 
centimeters, and V_a indicates the actual velocity in 
kilometers per hour. The machine and soil specific 
parameters as suggested in ASABE draft model are 
given in Table 2. 


Table 2 Machine Specific Parameters 


Implement Tillage Number of tines Machine parameter Soil Specific parameters 
A B C. Ff Fh Fs 
Primary 42 2.8 0 1 0.85 0.65 
Cultivator 7,9, 11... ete. 
Secondary 32 19 0 1 0.85 0.65 


2.4 Traction Performance Indices 

2.4.1 Motion Resistance Ratio (MRR) 
Motion resistance, also referred to as the motion 
resistance ratio (MRR), and represents the 
relationship between rolling resistance and the 
dynamic weight applied to the wheel [9]. 

2.4.2 Gross Traction Ratio (GTR) 
The gross traction ratio denotes the proportion of the 
tractor's gross thrust to the dynamic weight exerted 
on the traction wheels (powered wheels). 

2.4.3, Net Traction Ratio (NTR) 
The net traction ratio is the ratio of draft to the 
dynamic weight on the driving wheels. It can be 
computed using two methods, as outlined in 


equations 7 and 8. 


NTR, == 
7 


(7) 

NTR, = GTR— MRR 

2.4.4 Wheel Slip (%) 

When a tractor pulls a load, the distance traveled 
and/or speed decrease due to flexing of the tractive 
device and shear within the soil [10]. Slip occurs 
when a wheel or traction device experiences pull (net 
traction). It is given by equation 9. 

s=(1-%) x 100 (9) 

t 


Where s is slip (%), V,is actual velocity (km/hr) and 


(8) 
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V, is theoretical velocity (km/hr). 

2.4.5 Field Capacity (Ha/hr) 
It is defined as the actual area covered within a given 
time. 

3. Software Development in MATLAB Simulink 
Using MATLAB Simulink 2020a, a program was 
created to forecast slip, real operating speed, and field 
capacity. Figure 2 depicts the development process 
for the same Vis flow chart. The program was first 
fed with the input parameters. Input parameters were 
tractor static weight (Kg), weight of Implement (Kg), 
Hd (m), Xegt(m), Xegi(m), depth of operation(m), 
theoretical velocity, width of implement(m) and type 
of tire. The software first used equation 6 to 
determine the implement's draft based on the state of 
the soil. Next, every traction parameter was 
computed using equations 1, 2, 3, 7, 8, and 9. Motion 
resistance ratio (MRR) and slip were initially 
considered to be 0.04 and 0.02 respectively. They 
were then increased iteratively depending on if-else 
conditions. Based on the difference between MRR 
andMRR,., the MRR was increased and supplied as 
feedback to the program. Similarly, for slip 
difference between NTR, - NTR, was considered. 
After arriving at slip, the actual speed of operation 
and the field capacity can easily be calculated. 
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Figure 2 Flow Chart for Developed Software 
4. Instrumentation for Software Validation 
The tractor was equipped with standard tires sized 6- 
16 and 13.6-28 on the front and rear axles, 
respectively. Slip was determined using two 
measurements: the actual and theoretical speeds of 
the tractor. The theoretical speed was calculated 
based on the average rotational speed (rpm) of the 
rear wheels, while the actual speed was obtained from 
The front wheel. The front wheel had reluctor rings 
with 44 slots, while the rear wheel speed was 
measured by assessing the gear speed between the 
range gear and final drive. The rear reluctor ring was 
affixed to the rear axle shaft and secured to the flange 
using serrated wheel bolts, with the rear sensor 
bracket mounted at the rear axle center. Conversely, 
the front reluctor ring was mounted on the front 
wheel hub, and the sensor bracket was bolted onto the 
knuckle arm. Both configurations are depicted in 
Figure 3. Angular speed signals from the 
transmission gear and front wheel were generated 
through the transmission gear and reluctor ring at the 
front wheel, respectively. These signals, along with 
wheel speed sensors, were transmitted to a 
microcontroller for slip calculation. The signals 
received from the rear sensors were utilized to 
determine the theoretical speed by multiplying the 
average readings with the distance traveled by the 
rear wheel in one revolution on a hard surface. 
Similarly, the average output of the front wheel 
sensor was multiplied by the distance travelled by the 
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front wheel in one revolution to determine the actual 
forward speed. The wheel speed data were input to a 
microcontroller for slip calculation and digital 
display presentation. [11-15]. 


Pr 
Figure 3 Arrangement for Speed Sensor and 
Reluctor Ring 


1. Results and Discussion 

The developed model as shown in Figure 4 using the 
control algorithm as shown in Figure 2 had to be 
validated. To validate the same a 45 hp Mahindra 575 
DI tractor along with M&M Limited /360023 9 tine 
cultivator was used. The result of the experiment is 
shown in Table 3. Wherein, the operated gear was H1 
and M3 as recommended by manufacturer and the 
throttle was varied from 1600 to 1700. The cone 
index was 500+50 Kpa measured using cone 
penetrometer. The developed model predicted the 
slip within +4.03 to -4.55% relative variation [16]. 
((Predicted- experimental) x (100/Experimental)), 
actual speed of operation within +3.83 to -2.70 % 
relative variation and field capacity within +3.44 to - 
7.09% relative variation [17]. To predict the accuracy 
of the model, the predicted slip vs experimental slip, 
the predicted actual speed vs experimental actual 
speed, the predicted field capacity vs experimental 
filed capacity were plotted as shown in Figure 5 
((a),(b) and (c)). Results indicated a good correlation 
with coefficient of determination of (R*) 0.9727for 
slip, 0.9763 for actual speed and 0.9763 for field 
capacity. 
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Figure 4 Developed Model 


Table 3 Experimental Versus Predicted Data 


Slip (%) Actual Speed (Knv/hr) Field Capacity (Ha/hr) 
Gear | Throttle 
(rpm) Variation Variation Variation 
Experimental | Predicted | (%) Experimental | Predicted | (%) Experimental | Predicted | (%) 
© ® - ® @) a ® @®) 7 
E)*100E E)*100E E)*100E 

16.60 17.00 241 5.53 141 131 

Hi__| 1600 5.59 -1.07 -7.09 

3 si 2 2 

Hi__| 1700 ici ai . 6.03 ad -1.16 ae - 0.81 
17.00 17.10 059 6.08 127 1.24 

Hi__| 1800 6.23 -236 -2.36 

ga: ses 18.00 17.50 278 i 6.66 sae 133 136 fog 
19.00 18.50 2.63 6.76 138 139 

H1__| 2000 6.76 -0.07 0.72 

0.52 5 

eh ll aa 19.10 19.00 0.52 is 736 ane 145 1.50 ve 
19.20 19.00 -1.04 7.62 136 135 

Hi_| 2200 TAI 2.82 0.74 

eae 20.10 20.00 -0.50 eee 6.62 ree 143 141 58s 
21.00 20.50 -238 690 152 147 

M3__| 1700 7.00 “143 -3.29 
2 2 12 5 

Pa ee 1.00 20.50 38 ea 721 ax 158 1.56 Sone 
22.00 21.00 455 7.62 1.64 1.62 

M3__| 1900 12 -131 -1.22 
2 » 5 2 

ae 23.00 22.00 435 bas 8.00 ‘gay 1.75 172 pe 
22.00 21.00 455 8.40 1.78 1.76 

M3_| 2100 8.54 -1.64 “1.12 
2 2 3 

ae 24.00 23.20 333 7 8.60 ga 1.79 1.78 ase 
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Figure 5 Predicted vs Experimental Graph for Output Parameters 


Conclusion 

A model-based user-friendly software was developed 
for predicting slip, actual speed of operation and field 
capacity since they are predominant parameters to 
evaluate tractor-implement combination performance 
[18]. This study also includes the methodology to 
measure slip and validate the same by 
instrumentation which may be quite useful for 
researchers and engineers. The developed software 
was validated with experimental actual field results 
and found very close prediction for different 
operating parameters. 
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